Renormalization of Molecular Quasiparticle Levels at Metal-Molecule Interfaces: 

Trends Across Binding Regimes 
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When an electron or a hole is added into an orbital of an adsorbed molecule the substrate electrons 
will rearrange in order to screen the added charge. This results in a reduction of the electron 
addition/removal energies as compared to the free molecule case. In this work we use a simple 
model to illustrate the universal trends of this renormalization mechanism as a function of the 
microscopic key parameters. Insight of both fundamental and practical importance is obtained by 
comparing GW quasiparticle energies with Hartree-Fock and Kohn-Sham calculations. We identify 
two different polarization mechanisms: (i) polarization of the metal (image charge formation) and 
(ii) polarization of the molecule via charge transfer across the interface. The importance of (i) 
and (ii) is found to increase with the metal density of states at the Fermi level and metal-molecule 
coupling strength, respectively. 

PACS numbers: 85.65.+h,31.70.Dk,71.10.-w,73.20.-r 
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The position of an adsorbed molecule's frontier or- 
bitals with respect to the substrate Fermi level deter- 
mines the threshold energies at which electron trans- 
fer can take place across the metal-molecule interface. 
Such electron transfer processes represent a cornerstone 
of surface science and form the basis of photo- and non- 
adiabatic chemistry, organic- and molecular electronics, 
as well as scanning tunneling- and photoemission spec- 
troscopy [H, 0, H, [j, d, H H- Accurate descriptions of 
adsorbate energy spectra are thus fundamental for quan- 
titative modeling within these important areas. 

Recently, a number of experiments probing transport- 
and optical properties of molecules at metal sur- 
faces, have found strong reductions of electron addi- 
tion/removal energies due to polarization effects in the 
metal substrate 11 H, S S 11, see Fig. nil. Theoreti- 
cal studies of metal-molecule interfaces are usually based 
on density functional theory (DFT) and rely on an in- 
terpretation of Kohn-Sham (KS) eigenvalues as quasi- 
particle energies. This interpretation is unjustified in 
principle and questionable in practice - even with the 
exact exchange-correlation functional - because single- 
particle schemes like KS and Hartree-Fock (HF) take no 
account of dynamical effects such as screening of added 
electrons/holes. Inclusion of dynamical polarization ef- 
fects through classical image charge models have been 
applied to correct DFT calculations of quantum trans- 
port in molecular contacts Only recently, GW 
calculations by Neaton et al. have demonstrated that the 
energy gap of benzene is reduced by more than 3 eV when 
physisorbed on graphite due to image charge effects, jl^ 
As a complementary study, the present work focuses on 
qualitative trends and explores different binding situa- 
tions including strong-weak coupling and narrow-wide 
band substrates. 



To illustrate the problem, suppose we add an electron 
into the lowest unoccupied orbital (LUMO) of an ad- 
sorbed molecule. The associated energy cost is given by 
the quasiparticle (QP) spectral function 



where l^*^) and 1^'^''"^) denote the many-body A^- 
particle groundstate and {N + l)-particle excited states, 
respectively. The qualitative shape of yl^ is mainly dic- 
tated by the hybridization with the metal states. On the 
other hand, quantitative features such as the precise po- 
sition of its peak(s), depend on the system's response to 
the extra electron in the LUMO. This follows by noting 
that peaks in appear at dominant Fourier components 
of c\\^Q){t). Dynamical effects of this type can be in- 
corporated within an effective single-particle framework 
where the electron-electron interaction is represented by 
an energy-dependent potential known as a self-energy. 
In practice the self-energy must be approximated using 
e.g. many-body perturbation theory. For weak metal- 
molecule coupling, total energy calculations with con- 
strained occupations of the molecule represent a simple 
alternative to the many-body approach. 13J However, 
such calculations only provide a lower (upper) limit for 
the LUMO (HOMO) energy and, as we show here, these 
limits can differ from the true quasiparticle levels. 

In this Letter, we introduce a simple model of a metal- 
molecule interface and use it to investigate the effect 
of dynamical polarization on the molecular levels as a 
function of microscopic key quantities such as the metal- 
molecule coupling strength (F) and the metal density 
of states (DOS). We find that polarization effects con- 
sistently reduce the molecular gap as compared to the 
gas phase value. This effect is completely absent in HF 
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FIG. 1: (Color online) (a) Schematic of a molecule's HOMO 
and LUMO energy levels as it approaches a metal surface. 
For weak coupling (physisorbed molecule) the gap is reduced 
due to image charge formation in the metal. For strong 
coupling (chemisorbed molecule) dynamic charge transfer be- 
tween molecule and metal reduces the gap further, (b) The 
model used in the present study, (c) The semi-elliptical band 
at the terminal site of the TB chain and the resonances of the 
molecule. 



theory which systematically overestimates the gap. On 
the other hand, the KS spectrum is indirectly affected 
through the "exact density" requirement, however, the 
KS gap is always found to be too small. On the basis of 
GW calculations, we find that the gap reduction due to 
polarization of the metal, i.e. image charge formation, 
is insensitive to F but increases with the metal DOS at 
Ep. On the other hand, polarization of the molecule via 
charge transfer to/from the metal increases strongly with 
r leading to a direct correlation between adsorbate bond 
strength and the gap renormalization. 

Our model Hamiltonian consists of three parts, H = 
-ffmet + -f^moi + V , describing the metal, the molecule, 
and their mutual interaction, see Fig. [T|d. The metal is 
modeled by a non-interacting, semi-infinite tight-binding 
(TB) chain, 
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The molecule is modeled as two interacting levels repre- 
senting the highest occupied (HOMO) and lowest unoc- 
cupied (LUMO) molecular orbitals 

-ffmol = ^HflH + i£,H + ^ajflL + C^mol (3) 

U^oi = UofiH^nHi + UofiL^fiLi + UHLnHnL, (4) 



where e.g. fiH = c^h'\'^h\ + c^xCjyj, is the number oper- 
ator of the HOMO level. Notice that, despite the inter- 
actions, the eigenstates of Hmo\ are simply Slater deter- 
minants build from the orbitals \Ha) and \La). Finally, 
hybridization and interaction between the molecule and 
the terminal site of the chain is described by 
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E E *hyb(cj£,C^cr + cI^CQct) + UcxtSfloSN^ 



u=H,L cr=T,i 

(5) 

Here Sfio = {no — 1) and (5A^moi — {Nmoi — 2) represent 
the excess charge on the chain's terminal site and the 
molecule, respectively. We set Ep ~ corresponding to 
a half filled band, and adjust so that the molecule 
holds exactly two electrons in the groundstate. Specifi- 
cally this means — — Ao/2 — C/o/2 — Uhl- The model 
neglects interactions within the TB chain and between 
the molecule and interior TB sites (i < 0). These ap- 
proximations are, however, not expected to influence the 
qualitative trends described by the model. 

It is instructive first to consider the model in the 
limit thyb = (weak physisorption) where the metal 
and molecule interact only via the Coulomb term Ucxt- 
It is straightforward to verify that in this limit, the 
many-body eigenstates coincide with the HF solutions. 
More precisely they are single Slater determinants con- 
structed from the molecular orbitals, \Ha) and \La), 
and the eigenstates of the Hamiltonian, iJ,^^j(5A^moi) = 
-ffmot + Ucxt{SN„,o\)no, where (SNn^oi) = -2, . . . , 2 is the 
number of molecular orbitals in the Slater determinant. 
Although the eigenstates are single Slater determinants, 
the situation is different from the non-interacting case 
because the single-particle orbitals of the metal depend 
on the occupation of the molecule (each electron on the 
molecule will shift the potential of the terminal site by 
UcKt)- As we will see below this has important conse- 
quences for the position of the molecule's QP levels. 

In order to test the accuracy of the GW approximation 
for the problem at hand, we have compared it to exact 
diagonalization results in the case where the TB chain is 
truncated after the first two sites, see the supplementary 
material. [3] The GW spectral function is essentially 
identical to the exact result, also for strong interactions, 
demonstrating that GW captures the essential physics of 
the model accurately. 

To study how the molecular levels depend on the model 
parameters, we vary one parameter at a time keeping 
the remaining fixed at the following reference values: 
t = 10, ihyb = 0.4, Uo = 4, Uhl = 3, t/cxt = 2.8, 
Ao = 2. The reference values correspond to weak cou- 
pling (F « 0.02) and a wide-band metal (W ^ At = 40). 
The positions of the HOMO and LUMO levels are de- 
fined as the maximum of the corresponding spectral func- 
tions ^ which we calculate from the Green's function, 
An{e) = — (l/7r)ImG^„(£), see supplementary material 
for plots of a representative set of spectral functions 14 1 . 
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FIG. 2: (Color online) Position of the HOMO and LUMO 
levels as a function of different model parameters for a weakly 
coupled molecule (small F). Note that dynamical polarization 
is completely absent in IfF theory and incorrectly described 
by the "exact" KS theory. Open squares denote the exact 
total energy difference between the normal groundstate and 
the groundstate when the molecule has been constrained to 
hold one extra electron/hole. Lower right panel shows on- 
site elements of the static, linear response function of the 
HOMO, LUMO and terminal site of the TB chain. Clearly, 
the renormalization of the gap is due to polarization of the 
metal. 



The Green's function is obtained by solving the Dyson 
equation fully self-consistently in conjunction with the 
self-energy (I]gw[G] or Shf[G]) as described in Ref. [ist . 
Following the lattice version of DPT [l^ , we define the 
KS Hamiltonian as the non-interacting part of H with 
the on-site energies corrected to yield the exact occu- 
pation numbers. Using the GW occupations as "exact" 
target occupations this allows us to obtain the "exact" 
KS levels. 

In the upper left panel of Fig. [2] we show the posi- 
tion of the molecule's QP levels as function of the metal- 
molecule interaction, C/ext- The GW gap decreases as 
C/gxf The reduction of the gap corresponds to the en- 
ergy gained by letting the HP states of the metal relax 
in response to the added electron/hole, i.e. to the per- 
turbations ±[/ext^o- Indeed, the open squares show the 
difference in total energy between the normal ground- 
state and the groundstate when the molecule has been 
constrained to hold one extra electron/hole. [l7[ HP com- 
pletely misses this effect due to its neglect of orbital relax- 
ations, and consequently the HF gap becomes too large. 
In a screening picture, the difference between the HP and 



GW levels equals the binding energy between the added 
electron/hole and its image charge (the positive/negative 
induced density at the terminal site), corrected by the 
cost of forming the image charge. 

The "exact" KS theory underestimates the gap, but 
seems to reproduce the trend of the GW calculation. At 
first this is surprising since a mean-field theory cannot 
describe the dynamical effects responsible for the gap re- 
duction. The explanation is that the KS levels are af- 
fected indirectly: the KS levels are forced to follow the 
GW levels in order to reproduce the GW occupations. 

In the upper right panel of Fig. [H we show the depen- 
dence of the QP levels on the molecule's intrinsic gap, 
Aq. The fact that the renormalization of the QP gap 
due to image charge formation is independent of Ao fol- 
lows from the discussion above since HZ, is independent 
of Ao. 

It seems intuitively clear that the size of the gap reduc- 
tion (for fixed C/ext) should depend on the polarizability 
of the metal, i.e. how much the electron density changes 
in response to the perturbing field created by the added 
electron/hole. In the lower left panel of Fig. [2] we show 
the QP levels as function of the chain hopping, t.(l8| 
The gap reduction is larger for small t corresponding to 
a narrow band. This is easily understood by noting that 
a narrow band implies a large DOS at Ep which in turn 
implies a larger density response function. In the lower 
right panel of Pig. [5] we show the diagonal elements of the 
static component of the response function, XnnifjJ = 0), 
for the HOMO, LUMO and terminal site of the chain, 
respectively. In this weakly coupled regime, the response 
function of the HOMO and LUMO is clearly negligible 
for all values of t, while the response of the terminal 
site is significant and increases as t is reduced. This 
clearly demonstrates that the QP levels are renormal- 
ized by screening inside the metal. In general, only the 
response at frequencies |a;| < P is relevant as P^^ sets 
the decay time of the states c^lvl/o) and chI^^o)- 

Interestingly, the deviation between the GW levels and 
Ai?tot becomes significant for small t. We stress that this 
does not imply that the GW results are wrong. In fact, 
Ai?tot only represents an upper/lower bound for the true 
QP energies, and the deviation thus indicates that the 
overlap of c^lv^o) and cjy|^o) with higher lying excited 
states of the (A^-l-l)-particle system (see Eq. ([T])) is more 
important for small t. 

We now leave the regime of weak molecule-metal cou- 
pling and consider the dependence of the molecular spec- 
trum on ihyb, see Pig. [31 At finite thyb the molecular lev- 
els broaden into resonances, however, it is still possible to 
define the level position as the resonance maximum, (loj 
In addition to the GW, HP, and KS levels, we also show 
the result of a calculation where only C/moi is treated at 
the GW level while C/ext is treated at the HP level. This 
means that the only dynamical effects included in the 
GW(mol) calculation are those due to Umo\- In particu- 
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FIG. 3: (Color online) Left: Position of the molecule's HOMO 
and LUMO levels as a function of the hybridization strength, 
thyb- GW(mol) refers to a calculation where only the internal 
interactions on the molecule, f/moi, have been treated within 
GW while the interactions with the metal, (Xxt have been 
treated within HF. Right: Static response function for the 
HOMO, LUMO, and terminal site of the TB chain. The car- 
toon illustrates how dynamic charge transfer stabilizes the 
charged system and thereby reduces the gap. 



tween occupied and unoccupied molecular orbitals. The 
size of the gap reduction correlates directly with the 
static response function, and is promoted by larger metal 
and/or molecular DOS at the Fermi level. The latter sug- 
gests that transition metals with half-filled d band should 
be more effective in reducing the gap of molecular adsor- 
bates than e.g. alkali, sp, and the early/late transition 
metals. The results strongly indicate that dynamical po- 
larization is of fundamental importance for charge trans- 
fer at metal-molecule interfaces. 
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Lundbeck Foundation's Center for Atomic-scale Materi- 
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ETSF project. 



lar, image charge formation in the metal is ignored. 

Focusing on the GW and GW(mol) results, we con- 
clude that the reduction of the QP gap as function of 
thyb is mainly due to the interactions internally on the 
molecule. On the other hand, the reduction due to the 
molecule-metal interaction is largely insensitive to ihyb 
(this reduction is given by the difference between the GW 
and GW(mol) levels). From the right panel of Fig. [3l 
we see that \xnn\ for the HOMO and LUMO states in- 
creases with thyh indicating that the gap reduction due 
to Umoi is of a similar nature as the image charge effect, 
but with the molecule itself being polarized instead of 
the metal. Polarization of the molecule can occur via dy- 
namic charge transfer to/from the metal as illustrated in 
the cartoon of Fig.[3l Note that this picture is consistent 
with the fact that no polarization of the molecule occurs 
without coupling to the metal (see GW(mol) result in the 
limit thyb ^ in Fig. ^. 

The fact that dynamical polarization becomes more 
important as the molecular DOS at the Fermi level in- 
creases, has important consequences for charge transport 
in molecular junctions. Indeed, as the chemical poten- 
tial, fia, approaches a molecular level, polarization ef- 
fects become stronger and the level is shifted towards 
fia- This effect was recently shown to have a large im- 
pact on the junction IV characteristics. Finally, we 
mention that the correlation between chemisorption bond 
strength and gap renormalization suggested by Fig. [31 
has in fact been observed in inverse photoemission spec- 
troscopy experiments. 

In summary, we have used a simple model of a metal- 
molecule interface to identify universal trends in the way 
dynamical polarization renormalize molecular interface 
states. The effect of polarization is to reduce the gap be- 
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